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ABSTRACT 

We present 7 mm Mopra observations of the dense molecular gas towards the W28 
supernova remnant (SNR) field, following a previous 12 mm line survey of this region. 
These observations take advantage of the 7 mm beam size to probe the dense and 
disrupted gas in the region at ~1 arcmin scales. Our observations are focused towards 
the north-eastern (NE) HESS J1801-233 and southern HESS J1800-240B TeV gamma- 
ray sources, with slightly less observations towards HESS J1800-240 A & C. Using the 
CS (1-0) transition we reveal multiple regions of dense gas, nn 2 ~ 10 5 cm -3 . We report 
the discovery of dense gas towards HESS J1800-240C, at the site of a 1720 MHz OH 
maser. The NE molecular cloud is known to be disrupted, many 1720 MHz OH masers 
and broad CO line emission are detected at the rim of W28. Here we reveal this shock 
interaction region contains generally extended clumpy CS, as well as clumpy SiO and 
CH3OH emission with broad line profiles. The FWHM of the molecular lines extend 
up to 18kms _1 on the W28 side of the NE cloud. The detection of SiO towards maser 
clumps OH C, D, E & F provide further evidence of the shocked conditions in the NE 
cloud. Several other lines associated with star formation are also detected towards the 
southern source, notably the energetic Hll complex G5. 89-0. 39. The spatial match of 
dense gas with the TeV emission further supports the CR origin for the gamma-rays. 
We estimate the mass of several extended dense clouds within the field and predict 
the TeV flux from the dense cloud components. The predicted fluxes are on the order 
of 10 -14 -10 -13 phcm -2 s _1 , which should be detectable and possibly resolved by a 
future TeV instrument, such as the Cherenkov Telescope Array. 

Key words: ISM: clouds - Hll regions - ISM: supernova remnants 
- gamma-rays: observations - supernovae: individual: W28. 
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1 INTRODUCTION 

W28 is a striking example of TeV (10 12 eV) gamma- 
ray emission spatially ov erlapping with molecular 
gas (|Aharonian et all l2008bh . W28 is part of the in- 
creasing list of sources with spatial overlap between TeV 
or GeV (10 9 eV) gamma-ray emission a nd molecular gas 
(e.g. HESS J1745-290/SNR G35 9.1-0.5 (lAharonian et all 
2004), HESS J1714-385/CTB 37A (lAharonian et aZ.I|20Q8ah . 
HESS J1923+141/SNR G49.2-0.7 (iFeinstein et all I2009T L 
IC 443 ([Albert et a/.ll2QQ8l : lAcciari et aZ.ll2009h ), inaddition 
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to the central molecular zone ( CMZ) towards the Galactic 
centre region (|Aharonian et al\\200w . 

W2 8 is a mixed morphology old age (> 
10 4 yr; iKaspi et all [l993) supernova remnant (SNR) 
with dimensions of 50' x 45 / . W28 is e s timat e d to be at a 
distan c e of 1.2 to 3.3 kpc ( e.g. iGoudisI (|l976h : iLozinskaval 
(|l98lh : iMotogi et all (|201oh ). I t has been shown to exhibit 
non-thermal radio emission (iDubner et all l200d ). ther- 
mal X-ray emission dRho &; Borkowski I l20Q2h, and, more 
recent ly gamma-r ay emission at TeV (lAharonian et all 
l2008bh and GeV (|Abdo et all l20ld : iGiuliani et all l2010h 
energies. 

Several CO surveys in the (1-0), (2-1) and (3-2) 
lines reveal massive molecular clouds to the north-east 
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(NE) 
199a: 



and to the south (S) of the SNR (lArikawa et al, 



Reach et a l. 2005; T orres et a l. 2003; Aharonia n et al 



2008bl: iFukui et a/.ll2008r ). Most of the CO emission appears 
centred at a local standard of rest velocity (Vlsr) similar 
to that inferred for W28 Vlsr ~7kms _1 ( corresponding to 
a dis tance ~ 2kpc) based on Hi studies ([Velazquez et all 
l20ld ). 

The molecular clouds NE of W28 are known to ex- 
hibit broad CO line emission. It has been argued that W28 
has disrupted much of this gas, giving r i se to its broad ve- 
locity distribution (| Arikawa et al .1 ll999l : iReach et al. 2005; 
iTorres et all l2003h . The most likely mechanism for this 
type of disruption would be the SNR shock interacting 
with the molecular gas, as is the case f or other sources 
e.g. HESS J1745-290/SNR G35 9 1-0.5 (lAharonian et all 
l2004h . HESS J1714-385/CTB 37A (| Aharonian et alhOOSsl) . 
HESS J1923 +141/SNR G49.2-0 . 7 (iFeinstein et all Eooa ) 
and IC 443 (lAlbert et all 120081 : lAcciari et all 12009 ). All 
of these sources similarly display 1720 MHz OH masers 
and a re mature SNRs (age > 10 4 yr). Our NH3 observa- 
tions ([Nicholas et all l201Ch revealed that the W28 SNR 
shock has disrupted the dense core of the NE cloud, and 
strong NH3 (3,3) and (6,6) detections with broad line widths 
(FWHM > lOkrns -1 ), suggest that this region is warm and 
turbulent. The shocked NE cloud also provides an oppor- 
tunity to study the diffusion and propagation of CRs into 
molecular clouds. Additionally, the effects of a SNR shock 
propagating through a dense molecular cloud can be stud- 
ied and search for sites of star formation triggered by the 
passing shock. 

The southern clouds harbour sites of high mass 
star formation, c ontaining multi p le Hll regions, G6.225 - 
0.569, G6.1-0.6 ([Lockmanl Il989l : iKuchar fe Clark I Il997h 
an d the ultracompact Hll reg i on (UCHll) G5. 89-0 . 39, (see 
e.g lHarvev Forveillel (| 19881 ): iKim fe Kool (|200lh ). A re- 
cently detected 1720 MHz OH mase r towards a candidate 
SNR G5. 71-0.08 ([Brogan et all\200di ) may also suggest that 
there is another SNR shock with molecular gas interaction 
occurring in one of the southern clouds. A key question is 
whether the W28 SNR has disrupted the southern clouds 
and is responsible for the TeV emission, or whether other 
high energy processes are at play. The star formation activ- 
ity present in the southern clouds may play a part in the 
production of TeV emission, as theory into protost ellar par- 
ticle acceleration may suggest ([Araudo et all [20071 ). 

The the TeV gamma-ray emission detected towards 
the NE cloud has strong evidence supporting the hadronic 
emission mechanism. In order to model leptonic scenar- 
ios, uncomfortabl y low magnetic field s and gas densities 
are required ([Yamazaki et al | 2006l : iFuiita et all [2009; 
iGabici et all l20ld : lAbdo et all |201cF 7 These requirements 
contradict with the observed m olecular gas densities nu 2 > 
10 4 cm -3 ([Nicholas et a/.ll2010h and the m agnetic field en - 
hancements seen in dense molecular clouds (|Crutcherlll99lh . 
The key question relates to the origin of the TeV particles. 
Provided the hadronic emission scenario, a population of 
high energy particles are required to interact with the dense 
clouds. For the NE cloud, the SNR shock is the obvious 
source of accelerated particles. However, the same can not 
be conclusively said about the southern sources. CRs accel- 
erated by W28 in the past may have propagated south to 
illuminate the southern clouds. Additionally, multiple Hll 



regions and sites of star formation, including the energetic 
UCH11 complex G5. 89-0. 39, could assist in CR acceleration 
towards the southern TeV sources. 

Since the common and abundant gas tracer, CO, with a 
critical density of 10 2 cm -3 , rapidly becomes optically thick 
towards clumps and cores, probing the molecular cloud den- 
sity profile can be impaired. Ideal tracers of dense gas are 
those with a lower abundance (10 _5 xCO) and higher crit- 
ical density (10 4 -10 5 cm~ 3 ). Using the Mopra telescope, we 
(Nic holas et a/.ll201oh conducted broad scale (~ 1.5° square) 
observations of the W28 field in a 12 mm line survey. Taking 
advantage of the 23 GHz NH3 inversion lines, the cold dense 
interiors of the molecular clouds were probed and revealed 
dense gas spatially consistent with both the CO gas and TeV 
emission. 

In order to further probe the structure and details of 
the W28 field molecular clouds, we have continued to use the 
Mopra radio telescope in a 7 mm line survey to observe the 
CS (1-0) line, as well as other lines including SiO (1-0), Class 
I CH3OH masers and cyanopolyynes HC n N (n = 3, 5, 7), to 
trace the dense gas, the presence of shocks, outflows, and 
disrupted gas, and sites of high-mass star formation. 

The 43 GHz CS (1-0) line is a useful warm, dense gas 
tracer, and, as CS also exists in various isotopologue forms, 
the optical depth and column density can be constrained. 
The higher frequency of 7 mm lines compared to our previ- 
ous 12 mm observations, results in a smaller beam FWHM 
(rsj 1 arcminute). This makes CS an ideal follow-up tracer 
to determine the properties of the dense and disrupted gas 
toward the molecular clouds in the W28 region. 



2 MOPRA OBSERVATIONS AND DATA 
REDUCTION 

Observations were performed with the Mopra radio tele- 
scope in April of 2009 and March of 2010 and utilised the 
UNSW Mopra wide-band spectrometer (MOPS) in zoom 
mode. Mopra is a 22 m single-dish radio telescope located 
-450 km north-west of Sydney, Australia (31° 16' 04" S, 
149°05 / 59 // E, 866m a.s.l.). The 7mm receiver operates in 
the 30-50 GHz range and, when coupled with MOPS, allows 
an instantaneous 8 GHz bandwidth. This gives Mopra the 
ability to cover 40% of the 7 mm band and simultaneously 
observe many spectral lines. The zoom mode of MOPS al- 
lows observations in up to 16 windows simultaneously, where 
each window is 137.5 MHz wide and contains 4096 chan- 
nels in each of two polarisations. At 7 mm this gives MOPS 
an effective bandwidth of — lOOOkms -1 with resolution of 
—0.2 km s _1 . Across the whole 7 mm band, t he beam FWHM 
varies from 1.37' (31 GHz) to 0.99' (49 GHz) ([Urauhart et all 
l2010h . Table □ lists the lines which MOPS was tuned to re- 
ceive. 

'On-the-fly' (OTF) mapping observations were based 
on our earlier 12 mm mapping of the W28 region. Although 
we mapped all four TeV sources, deeper observations to- 
wards the NE shocked cloud (HESS J1801-233) and south- 
ern cloud (HESS J1800-240 B) were obtained than toward 
HESS J1800-240 A & C. The observations were recorded al- 
ternating the scanning direction from RA to Dec. to reduce 
noise levels and to eliminate artificial stripes that can be 
introduced when only one scanning direction is used. Col- 
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Table 1. Molecular lines and their corresponding rest frequencies 
which MOPS was tuned to receive. The final two columns indicate 
whether the line was detected in our mapping or deep pointing 
observations. 



Molecular Line 


Frequency 


Detected 


Detected 


Name 


(MHz) 


Map 


Deep Spectra 


— 

30 SiO (1-0, v = 0) 


42373.365 


~ 




SiO (1-0, v = 3) 


42519.373 


Yes 




SiO (1-0, v = 2) 


42820.582 


Yes 


~~ 


29 SiO(l-0, v = 0) 


42879.922 




-~ 


SiO (1-0, v = 1) 


43122.079 


Yes 


— 


SiO (1-0, v = 0) 


43423.864 


Yes 


Yes 


CH3OH-I 


44069.476 


Yes 


Yes 


HC7N (40-39) 


45119.064 






HC 5 N (17-16) 


45264.75 


Yes 


Yes 


HC 3 N (5-4) 


45488.839 


Yes 


Yes 


13 CS (1-0) 


46247.58 


Yes 


Yes 


HC 5 N (16-15) 


47927.275 


Yes 


Yes 


C 34 S (1-0) 


48206.946 


Yes 


Yes 


OCS (4-3) 


48651.6043 






CS (1-0) 


48990.957 


Yes 


Yes 



lectively, over both observation runs we obtained 9 passes on 
the NE cloud (HESS J1801-233), 9 passes towards the south- 
ern cloud (HESS J 1800-240 B) and 3 passes toward each of 
the HESS J1800-240 A and HESS J1800-240 C regions. 

Data were reduced using the standard ATNF packages 
livedata, gridzilla, ASAP and miriadQ For mapping data, 
LIVEDATA was used to perform a bandpass calibration for 
each row using the preceding off-scan as a reference, and 
applied a linear fit to the baseline, gridzilla re-gridded and 
combined all data from all mapping scans into a single data 
cube, with pixels (Ax, Ay, Az) = (15", 15", 0.21 kms" 1 ). The 
mapping data were also weighted according to the relevant 
Tsys, Gaussian-smoothed based on the Mopra beam (1.2' 
FWHM and 3' cut-off radius), and pixel masked to remove 
noisy edge pixels. Data cubes were then converted into an 
ASAP scantable and had a 7 th order polynomial function 
fit subtracted from each pixel in the cube before returning 
the final output fits files. Analysis of position-switched deep 
pointings employed ASAP with time- averaging, weighting by 
the relevant Tsys, and baseline subtraction using a linear 
fit after masking the 15 channels at each bandpass edge. 

In both mapping and position-switched data, the an- 
tenna temperature T^ (corrected for atmospheric attenu- 
ation and rearward loss) was converted to the main beam 
brightness temperature T m b (K), such that T m b = T\/r] m b 
where r] m h is the Mopra main beam efficiency. Eac h line 
was co rrected by the relevant f] m b following lUrquhart et all 
(|2010h . 

Due to the uneven mapping exposure across the W28 
field, the RMS error in T m b, Trms (K), varies. On top of the 
exposure dependence, there is also a frequency dependence 
on Trms across the 42-49 GHz band. For this reason in Fig- 
ure [1] we present a CS (1-0) peak pixel map of the W28 field 
with dashed boxes outlining the exposure boundaries. The 



1 See |http : / /www . atnf . csiro . au/ computing/ software/ for 

more information on these data reduction packages. 



numbers presented in the white boxes are the averaged Trms 
achieved (across the 42-49 GHz band) within that boundary. 
We also show in Figure [2] the velocity of the peak pixels to 
show the general velocity locations and distribution of the 
various clumps of CS (1-0) gas. We show the same contours 
on both Figures [1] and [2] for guidance. 

For the detected lines, images of velocity-integrated 
intensity, position velocity (PV) and velocity dispersion 
(vrms) have been produced where possible (although not 
all are shown). Velocity-integrated intensity images (miriad 
th moment) are integrated over a velocity range determined 
by the observed line width. This was chosen to encompass 
the bulk of the emission from the region of interest. Mini- 
mum contour levels were set based on the integrated Trms 
of the emission. The integrated Trms of each image was de- 
termined by creating additional MIRIAD th moment maps in 
a velocity space either side of the velocity range of interest, 
ensuring the same number of channels were used. These ad- 
ditional moment maps were used to create pixel distribution 
histograms which were Gaussian fitted. In this way we es- 
timated the integrated Trms from the Gaussian fit results. 
On all images, the minimum contour levels are mentioned in 
terms of the raw value as well as Trms- Generally the min- 
imum accepted contour level on integrated maps is 2 Trms • 
PV plots were created by re-ordering the data cubes axes 
and Hanning smoothing the velocity axis (width ~ 3 km s _1 ) 
to improve image quality. The PV plots show the peak pixel 
along the declination axis for illustrative purposes. Inten- 
sity weighted velocity dispersion (miriad 2 moment), vrms, 
maps were calculated for pixels above a reasonable thresh- 
old. 



3 ANALYSIS AND RESULTS OVERVIEW 

Of the seventeen molecular transitions searched for (Ta- 
ble [1]), eleven were detected. From the mapping data we 
detected: CS, C 34 S and 13 CS in the (1-0) transition; SiO (1- 
0) in several vibrational modes; CH3OH; HC3N (5-4) and 
HC 5 N in the (17-16) and (16-15) transitions. 



3.1 CS Emission 

The CS (1-0) line was detected across the W28 field and 
is the most prominent. Along with the common CS iso- 
topologue, 12 C 32 S we also detect other isotopologues, C 34 S 
and 13 CS in the (1-0) transition. Unless otherwise indi- 
cated, CS refers to the common 12 C 32 S isotopologue. From 
our deep pointing observations (DP-1 to DP-5, locations 
shown in Figure [T]) with Trms 00 2 times lower than the 
mapping data, a Gaussian fit to the detected isotopologue 
line emission allowed the optical depth and column den- 
sity to be estimated. The method assumed an elemental 
abundance ratio between CS and the rarer isotopologues 
of 22.5 for [CS]/[C 34 S ] and 75 for [CS]/[ 13 CS 1 to calcu- 
late the optical depth using Zinchenk o et al. I (|l994l . Equa- 
tion 1). The upper st ate column density wa s then estimated 
via the method from iGoldsmith fe Langerl (l999, Equation 
9). The total column density of all states was obtained by 
applying a local thermal equilibrium (LTE) approximation 
and usin g temperature estim ations from our NH3 obser- 
vations ([Nicholas et q/.ll201Ch . assuming that the CS gas 
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Figure 1. CS (1-0) peak pixel map of the W28 region (log scale image with white linear contours, min. contour ~ IOTrms) with 
black/white crosses + to indicate the locations of the position-switched deep pointings (DP-1 to DP-5). Thick black dashed contours are 
the H.E.S.S. TeV emission (4a, 5cr, 6cr levels) revealing HESS J1801-233 and the different components of HESS J1800-240. Boundaries 
of catalogued SNRs are indicated by dashed circles, and locations of Hll regions are indicated by black/white crosses x. Black dashed 
boxes are used to separate regions based on their exposure. The numbers presented in the white boxes are the mean main beam Trms 
(units of K) achieved across the 42-49 GHz band within that boundary. The 1 arcmin beam FWHM is also shown in the lower left corner. 



is tracing the same gas as our NH 3 observations. From 
these parameters, the LTE hydrogen mass, number density 
a nd virial mass were e stimated via the methods outlined 
in lNicholas et al . (2010), by assuming a compact-source size 
r = 0.2 pc and a molecular abundance ratio, Xcs, of CS to 
H2. CS molecular abundance ratios are known to vary dra- 
matically f rom different s i tes in the Galaxy, ranging from 
1 x 10" 8 (iFrerking et all Il980l, eg . Sgr B - Table 4) to 
6 x 10" 11 (|Linke fc Goldsmith II 1980 . eg. W43 - Table 3). We 
adop ted an abu ndance ratio of Yes — 4 x 10 _ 9 (|lrvine et al\ 
Il987h following IZinchenko et aU (|l994h which is typical for 
dense quiescent gas. This is likely the case for most of the 
clumps, with the exception of the NE cloud as later dis- 
cussed in section 15.11 The results of our compact-source size 



LTE analysis and cloud mass/densities from deep pointing 
observations are displayed in Table [2] 



The LTE analysis method has also been applied to 
larger extended regions by utilising spectra from mapping 
data. For extended source calculations, we extracted and 
averaged the spectra for all pixels contained within an el- 
liptical region. In this way, we apply the standard analysis 
outlined above to a single averaged spectra for the entire 
region. We used the same abundance ratio, but assumed the 
emission fills a volume denned by a prolate ellipsoid of vol- 
ume V — 4/3 x 7rrir%. Results of the extended source LTE 
analysis from mapping data are displayed in Table [3] 
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Figure 2. Velocity of the peak pixels, with the same black and white contours as shown in Figure ^ Red solid ellipses are the regions 
where spectra were averaged in order to extract gas parameters as an extended source. 



3.2 SiO Emission 

We also tuned MOPS to include several different SiO (1- 
0) transitions arising from different vibrational modes (v = 
0,1,2,3), as well as the isotopologues 29 SiO and 30 SiO. 
The SiO molecule is a known and prov en signpost for the 
presence of shocked and disrupted ga s (Schilk e et al. 1997; 
iMartin-Pintado et aZ j2005lGusdorf et a l. 2008a, b). Our de- 
tection of SiO from several locations provides some insight 
into the disruptive conditions within the molecular clouds. 
Unless otherwise indicated, SiO refers to 28 SiO in the v = 
mode. 



3.3 CH 3 OH Emission 

Multiple detections of CH3OH masers towards the southern 
clouds are indicative of the star formation activity present 
towards the Hn regions and surrounding environments. Ad- 
ditional detections of CH3OH emission in the NE cloud is 



interesting. To date, no star formation tracers such as: H2O 
masers, CH3OH masers or infra-red emission have been de- 
tected towards this cloud. The CH3OH emission in the NE 
cloud does not appear to correspond to stellar activity, as 
discussed later. 



3.4 Cyanopolyyne HC n N (n = 3, 5, 7) Emission 

Detections of the cyanopolyynes HC 3 N (5-4), HC 5 N (16-15) 
and HC5N (17-16) towards the southern clouds are indica- 
tive of activity in cold dense cores at the very early 
stages of high mass star formation. Interestingly, ad- 
ditional detections of HC3N (5-4) towards the shocked NE 
cloud are revealed. The sites of HC3N (5-4) emission in the 
NE cloud do not appear to correspond to stellar activity. 

A more detailed discussion of the morphology and im- 
plications of the molecular transitions detected from each of 
the mapped regions follows in Section [H 



© 0000 RAS, MNRAS 000, 000-000 



6 B. Nicholas et al 



Table 2. Calculated gas parameters from CS (1-0), C 34 S (1-0) and 13 CS (1-0) isotopologue ratios from deep pointing observations DP-1 
to DP-5 (Figured}. Columns from left to right are: region name, CS (1-0) peak T mb temperature, CS (1-0) LSR velocity, CS (1-0) FWHM, 
CS (1-0) optical depth r, assumed kinetic temperature T^, beam corrected total column density iV[CS], hydrogen column density AffH^], 
mass, molecular hydrogen number density tih 2 an d virial mass range M w - lY . The virial mass range is the lower and upper bounds when 
considering a constant, Gaussian and r~ 2 density profile. 



Core / 


T m b 


Vlsr 


FWHM 


T 


T k * 


7V[CS] 


N[R 2 ] 


Mass 


n H 2 


Virial Mass 


Region 


[K] 


[kms- 1 ] 


[kms" 1 ] 




[K] 


[xl0 14 cm- 2 ] 


[xlO 23 cm" 2 ] 


[M ] 


[xl0 5 cm- 3 ] 


[xlO 2 M ] 












- Compact Source - 










DP-1 


1.1 


13.1 


18.1 


1.0 


46 


16.5 


4.1 


790 


6.9 


133 - 279 


DP-2 


2.4 


7.4 


11.2 


1.6 


46 


28.8 


7.2 


1400 


12.0 


23 - 81 


DP-3+ 


1.4 


9.0 


2.2 




20 


0.1 


2.2xl0~ 2 


77 


0.5 


1.2 - 4.3 


DP-4 


9.5 


9.3 


4.0 


2.8 


33 


61.1 


15.3 


2900 


25.4 


3.4 - 11.9 


DP-5* 


2.8 


13.0 


1.1 


2.1 


20 


4.0 


1.0 


190 


1.7 


0.3 - 0.9 



Compact-source mass/density scaling factors vs. radius R (pc) 



R (pc) 

Mass M(R)/M (0.2 pc) 
Density n H2 (R)/n H2 (0.2 pc) 



0.10 0.15 0.20 0.25 0.30 0.35 0.40 
0.89 0.93 1.00 1.09 0.21 1.36 1.54 
7.09 2.21 1.00 0.56 0.36 0.25 0.19 



* Temperatures from NH3 observations (Nicholas et aE| |201Cl|) . 
t Using only CS (1-0) and assuming optically thin emission. 

* Using only CS (1-0) and C 34 S (1-0). 



Table 3. Calculated gas parameters from CS (1-0), C 34 S (1-0) and 13 CS (1-0) isotopologue ratios for extended sources, using averaged 
spectra from mapping observations. Columns from left to right are: region name, peak T mb temperature, CS (1-0) line FWHM, optical 
depth r, assumed kinetic temperature Tk, beam corrected total column density iV[CS], hydrogen column density iV[H2], mass and 
molecular hydrogen number density tih 2 • We treat the NE shocked cloud mass as an upper limit resulting from uncertainties in the 
molecular abundance ratio. 
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a For ellipse 7.2 x 2.8 pc diam.; 


pos. an 


gle 


-30°; centred on RA 18:01:46.7 Dec -23:24:21.9 


shown in Figure [2] 






b For ellipse 4.8 X 4.2 pc diam.; 


pos. angle 


+ 30°; centred on RA 18:00:32.9 Dec -24:03:59.2 shown in Figure[2] 






c For ellipse 6.3 X 2.1 pc diam.; 


pos. an 


gle 


-50°; centred on RA 18:00:22.0 Dec -23:56:59.2 shown in Figure[2] 






d For ellipse 7.0 X 2.8 pc diam.; 


pos. an 


gle 


-30°; centred on RA 18:00:55.9 Dec -24:09:14.1 


shown in Figure [2] 






e For ellipse 5.6 X 2.8 pc diam.; 


pos. an 


gle 


-10°; centred on RA 18:01:45.2 Dec -23:59:51.2 


shown in Figure [2] 






f For ellipse 2.8 X 1.4 pc diam. 


pos. an 


gle 


+ 30°; centred on RA 18:02:00.7 Dec -23:53:01.4 


shown in Figure [2] 







4 SPECIFIC DISCUSSION OF REGIONS 

4.1 North-eastern Cloud / HESS J1801-233 

Located at the northern boundary of HESS J1801-233 and 
nearby the giant Hn region M20, we observe several molecu- 
lar lines, which reveal dense gas and star formation activity. 
These molecular lines are all at a Vlsr range ~ 20 kms -1 , 
which is further away than the rest of the molecular emis- 
sion seen in the NE cloud, Vlsr ~ 10 kms -1 , discussed later. 
This core was detected in our 12 mm study and where it was 
referred to as Core 1 . Toward t his northern most core is an 
IRAS source, IRAS 17589-2312 (iBronfman et aZj|l996h , and 
we adopt the IRAS nomenclature. The dense CS core seen 
here is likely part of a molecular bridge connecting down to- 
wards the peak of the HESS J1801-233 TeV emission seen in 



Figure [T] We note that in the peak pixel map, Figure [T] the 
lowest contour level is ~ IOTrms- This gives an indication 
to the prominence and the extent of the CS (1-0) emission 
detected towards the NE cloud and the W28 clouds in gen- 
eral. Toward IRAS 17589-2312, we observe multiple Class 
I CH3OH masers along the line of sight within the beam. 
The maser emission is relatively weak (peak T mD ^0.6 K) 
compared to the other CH3OH masers seen in the W28 field 
discussed later. This detection adds to the list of maser tran- 
sitions w hich have been obse rved from this r egion: 6.7 GHz 
CH3OH (iFontani et q/Jl201Ch . 22 GHz H 2 (ICodella et all 
ll995l : lNicholas et q/Jl201oF aiid now 44 GHz CH 3 OH. 

We also detect the cyanopolyyne HC3N (5-4) line to- 
wards this core indicating hot gas phase chemistry probably 
driven by the central source. There is also a detection of 
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SiO (1-0) which indicates s hocks and /or outflow s originat- 
ing from IRAS 17589-2312. iLefloch et all <|2008h label this 
core as TC5 and suggest that the star formation in this core 
is recent and could possibly have been triggered by the W28 
SNR. Spectra for the line emission detected towards this 
core are shown in Figure [ATI 



Further south an extended dense region is spatially con- 
sistent with the TeV gamma-ray source HESS J 180 1-233. 
Towards the north-east rim of W28 several mm-line sur- 



vevs (IWoottenlll9~8ll: 


Arikawa et aZ.lll999b Reach et aZ.ll2005l: 


Torres et al 120031: 


Aharonian et all l2008bt Fukui et all 


20081: Nicholas et al. 


|2010h have revealed the well-known 



SNR shock interaction site, traced by broad line mea- 
surement s and the presence of many 1 720 MHz OH 
masers ([Frail et all\ 19941 ; Iciaussen et aZ.lll997h . 

Figure [3] contains integrated intensity images for emis- 
sion detected towards this region. The CS (1-0) emission is 
spatially matched to t he CO (2-1) emis sion measured by 
the Nanten telescope (jFukui et q/J 12008). and NH3 emis- 
sion measured by Mqpra r evealed a dense clump in the NE 
cloud ([Nicholas et q/jEoiol ) . This spatial overlap is generally 
expected as CS will trace the denser (n ^10 5 cm -3 ) regions 
of the molecular clouds. An image showing the spatial rela- 
tionship between the CO (2-1), CS (1-0) and NH 3 (3,3) gas 
is the NE cloud is shown in Figure IA3I 

The arcmin beam size allows the dense portion of the 
cloud to be probed in more detail than in our two arcmin 
beam NH3 study. Figure [3] highlights the clumpy nature of 
the cloud. Several dense clumps are arranged in a north-east 
to south-west direction. We also note the minimum contour 
level for the integrated CS (1-0) emission in Figure [3] is set to 
IOTrms- The CS (1-0) emission is very prominent through- 
out the NE cloud and likely extends towards the limits of 
our mapping. 

The line profiles of the CS and isotopologue emission 
are broad, with FWHM ^18kms -1 consistent with other 
molecular line observations from the region. Additionally, 
the asymmetry in the spectra provi de evidence for th e pres- 
ence of shocks and gas disruption ([Seta et aZ.I Il998l ) . Some 
CS line profiles indicate wings on the 'positive', or far side of 
the peak emission (e.g. CS (1-0) point 3/DP-l on Figure [3} 
and also on the 'negative', or nearer, side of the peak emis- 
sion in others (e.g. CS (1-0) point l/DP-2 on Figure [3}. In 
addition to the asymmetric line profiles, symmetric line pro- 
files (or at least profiles with similar wings on both sides of 
the peak emission) are also seen towards this NE cloud (e.g. 
CS (1-0) point 2 on Figure [3}. These line profiles may sug- 
gest there are components of the shock travelling in different 
directions, the possibility of reverse shocks or may simply be 
a projection effect. 

Several clumps of SiO (1-0) are detected and bounded 
by an arc of 1720 MHz OH masers, further confirming that 
a shock has disrupted this molecular cloud. The peak of 
the SiO (1-0) emission (point 15 in Figure [3} lies offset from 
the peak of the dense gas traced by CS (1-0) (point 2 on 
Figure [3}, which may suggest that the densest CS (1-0) core 
is shielded from the disruption of the shock. 

A class I CH3OH maser transition is detected toward 
the eastern edge of the dense gas (point 13 on Figure [3}. 
This is the weakest of the CH3OH masers detected, with a 
peak temperature of ~ 0.4 K. A second possible site of emis- 
sion appears real from the DP-2 spectra (point 14/DP-2 on 



Figure [3}. Here the spectra is weak, peak T m b 0.06 K, but 
with broad FWHM 15kms _1 . Both these CH3OH emission 
regions are not associated with typical star forming cores. 
Rather, they appear to be a by-product of the SNR shock 
heating the gas in a similar manner as di scussed, for the 
case o f expanding shocks from Hli regions in lVoronkov et all 
(|2010h . Similarly, the detection of HC 3 N(5-4) from this re- 
gion further indicates hot core gas phase chemistry. The 
emission is again broad, FWHM -12 km s" 1 as for other 
lines detected from the region. 

The position velocity (PV) plot shown in Figure H] high- 
lights the CS (1-0) emission, especially from the western or 
W28 side of the dense cloud. The velocity at the peak pixel 
map shown in Figure [2] reveals the existence of a velocity gra- 
dient, where the peak emission towards the south-western 
edge of the cloud has larger Vlsr values which decrease 
towards the north-east. This gradient would suggest that 
the shock is disrupting the cloud from the south-western or 
W28 side. The velocity dispersion vrms image for the CS (1- 
0) emission (for T m b >0.5K in the -5 to 20kms _1 range) 
from the north-eastern region is also shown in Figure 2) The 
emission with the largest v-rms dispersion originates from 
the south-western side of t he cloud, which was a lso the case 
for NH3 (3,3) emission in iNicholas efall (|2010h . The more 
disrupted emission is generally offset from the peaks seen in 
the integrated map, which are also shown on the PV plot 
(Figure [4} as white contours. Interestingly one of the clumps 
seen in the CS (1-0) integrated map (Point l/DP-2 in Fig- 
ure [3} appears disrupted and has a large v-rms value. This 
clump also corresponds to a peak in the SiO (1-0) integrated 
emission (Point 17/DP-2 in Figure [3}. 

We have taken two position switched deep pointings to- 
wards the NE cloud, specifically targeted towards the west- 
ern side of the cloud where the most turbulent conditions 
are seen. The locations of the position switched observa- 
tions are illustrated in Figure [1] DP-1 is towards the broad- 
est and most disruptive region, and DP-2 is towards a dis- 
rupted core. From these deep spectra, we have been able 
to extract gas parameters, which are displayed in Table [2] 
We find that the compact-source size mass, assuming a core 
with radius r = 0.2 pc, is more than a factor 10 smaller than 
the mass obtained when using the virial theorem. This is to 
be expected as these observations are targeted to the most 
disrupted regions of the cloud which are not expected to 
be gravitationally bound. Additionally, the molecular abun- 
dance ratio may be quite different in this shocked region 
compared to the unshocked gas. 



4.2 South-eastern Cloud / HESS J1800-240A 

The region associated with HESS J 1800- 24 A contains two 
Hn regions, G6 .225-0 .569 and G6.1-0.6 ([Lockmanl [l989l : 
iKuchar fc Clark Ill997h and IRAS 17588-2358. Integrated in- 
tensity maps of detected emission are shown in Figure [5] 
Dense gas traced by CS (1-0) populates this region in clumpy 
cores near the two Hn regions. A core unrelated to the Hn 
regions (point 3 in Figure [5} shows signs of extension and 
also displays a double peaked line profile towards its edge 
(point 4 in Figure [5}. Further south, additional clumps are 
detected with > 5 integrated Trms- The line profiles of the 
CS emission towards these cores are also moderately broad, 
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Figure 3. Integrated intensity images and spectra for the north-eastern HESS J 1801-233 region. Integration ranges are indicated on 
each panel. In all panels, H . E.S.S. TeV emission is indicated by thick black contours (4 &; 5cr levels), the positions of 1720MHz OH 
masers from Clauss en et al ] (I1997T ) are indicated by blue crosses +, and white contours are used to highlight the molecular emission 
from the image. Minimum contour level for: CS (1-0) is 2.6Kkm S - 1 (~10T RM s), for C 34 S (1-0) is OiKkms" 1 (~2T RMS ), for 13 CS 
(1-0) is 0.24Kkms- 1 (-2T RMS ), for SiO (1-0) is 0.32Kkm S - 1 (~2T RM s), for CH3OH ClKkms" 1 (~2T RMS ) and for HC 3 N(5-4) 
is 0.5Kkms — 1 (~2 Trms)- Surrounding the integrated intensity images are spectra taken from the both the mapping data and deep 
pointings, labelled according to their locations. 
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Right Ascension (J2000) Right Ascension (J2000) 



Figure 4. Left: Position velocity (PV) plot for the north-eastern region CS (1-0) emission. This RA vs Vlsr image is showing the peak 
pixel along the declination axis. Black /white crosses + represent the positions of the 1720 MHz OH masers with their corresponding 
labels from Claussen et al 1 (I1997T ). The separate core at Vlsr ~ 20kms 1 corresponds to the IRAS source. Right: Velocity dispersion 
urms map for the CS (1-0) emission from the north-eastern region. The ^rms nas been calculated for emission ^ 0.5 K in the -5 to 
20kms _1 interval and then scaled to show the line FWHM. White contours are the integrated CS contours seen in Figure [3] In both 
panels the RA axis is shown in degrees for clarity. 



FWHM rsj 5kms _1 , but not to the extent of the emission 
seen in the NE cloud. 

In the NE corner of our mapping, near the Hn region 
6.225-0.569, we trace star formation activity, detecting a 
bright Class I CH 3 OH maser (peak T mb -18 K), a HC 3 N (5- 
4) emitting core and evidence of shocks/outflows with a de- 
tection of SiO(l-O). The location of these star formation 
tracers are offset from the catalogued position of G6.225- 
0.569. 

We note that this region, associated with HESS J1800- 
240 A, has less exposure (a factor 3 lower than the north- 
eastern HESS J1801-233 and southern HESS J1800-240B 
regions) . 



4.3 Southern Cloud / HESS J1800-240B / 
G5.89-0.39 

Spatially well matched to HESS J1800-240 B are three dense 
cores oriented in a south-east (SE) to north-west (NW) 
alignment (points 1, 2 & 3 in Figure HQ), which were seen 
in our previous 12 mm observations ([Nicholas et all l2010l 
labelled as Triple Core SE, Central and NW). The SE and 
central sources are associated with the IR-bright and ener- 
getic UC Hn region G5. 89-0. 39 and the NW source with a 
pulsating M type star, V5561 Sgr, a further 5 arcmin away. 

The extraordinary molecular outflow emanating from 
the Hn complex G5. 89-0. 39 has been extensively stud- 
ied from a rcsec to arcmin sca l es in m ultiple molecular 
lines (e. g. lHarvev &; Forveillel (Il988h ; IChurchwell et al 



1990h : IZiilstra et al 



lAcord et al 



Gomez f (fl99lh: IChoi e~l 



Claussen et al 



Thompson fe Macdonaldl (Il999h: iKim fe Kool feOOlLkoodi: 
Sollins et all (|2004h : iHunter et all (|2008T ): iNicholas et all 
42Q10|)) 

Outflow velocities ±70kms -1 are seen in CO lines and 



similar claim s are m ade for SiO lines by lAcord et al ■1 (|l99/t ). 
Sollins et al argues that the outflow originates from 

a 1.3 mm continuum s ource and not from Feldt's candidate 
star (jFeldt et all 12003), as the star is not equidistant from 
the outflow lobes. The Hll complex G5. 89-0. 39 is actually 
comprised of two active star formation regions, Hn G5.89- 
0.39 A to the east, and UCHll G5.89-0.39 B approximately 
2 arcmin to the west (|Kim fe Kooll200lh . G5.89-0.39A is 
also known as W28A2 with strong radio continuum emis- 
sion, while G5.89-0.39B possesses the outflow. We observe 
SiO (1-0) over a wide region towards the three sources sug- 
gesting there are shocks or outflows. Multiple higher order 
SiO transitions on smaller scales are detected towards the 
G5. 89 -0. 39 B region and have been discussed by lAcord et all 
(E993). 

A new perspective that this work provides is the broad 
scale (15 arcmin) dense gas surrounding the Hn regions. NW 
of the UC Hn region G5.89-0.39 B in CS (1-0) we see an ex- 
tended arm feature with length ~ 6 arcmin extending north 
towards W28 (e.g points 3 & 4 in Figure [6}. The dense gas 
seen in this arm appears to have a slight velocity gradient, 
approximately +0.3kms _1 arcmin -1 from north to south. 
To the SE of the Hn regions there are several dense clumps 
which trace the dense com ponent of the cloud seen in CO (2- 
1) bv iFukui et~al\ <|2008h . Once again, these cores display 
broad CS (1-0) emission (FWHM>4kms _1 ) however, these 
clumps are at a slightly larger Vlsr compared to the Hn re- 
gions. Interestingly, most of the CS (1-0) emission from the 
southern cloud contains multiple clumps along the line of 
sight, all within ^20kms -1 of each other. This makes i>rms 
calculations difficult as the multiply clumped emission pro- 
duces artificially broad results. 

Star formation tracers such as 44 GHz class I CH3OH 
masers and HC3N (5-4) emission have been detected towards 
all three of the aligned sources. We also detect signs of warm 
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Figure 5. Integrated intensity images and spectra for the emission seen towards HESS J1800-240A. Integration ranges are indicated 
on each panel. In all panels, black contours are the H.E.S.S. TeV emission (4-6<j), black/white crosses x indicate the positions of the 
Hll regions G6. 225-0. 569 (north) and G6. 1-0.6 (south), the cross + is the location of IRAS 17588-2358, and the white contours are used 
to highlight the molecular emission from the image. Minimum white contour level for CS (1-0) is 1.75 K km s -1 (~5 Trms); f° r C 34 S 
is 0.52Kkms- 1 (~2T RMS ); for HC 3 N (5-4) is l^Kkms" 1 (~4T RM s); for CH3OH is 0.22Kkm S - 1 (~3T RM s) and for SiO (1-0) is 
0.42Kkms _1 (~3T R ms)- Surrounding the integrated images are spectra taken from the mapping data and are labelled according to 
their locations. 



gas phase chemistry via detections of HC3N (5-4) towards 
dense cores SE of the Hn regions (e.g point 22 in Figure [6]) 
which lie in a dense molecular ridge (see Figure [T]). We also 
detect the HC 5 N (16-15) and (17-16) lines towards G5.89- 
0.39 B. This would indicate that there is considerable energy 
available to excite the molecules into such a highly excited 
state. 

An additional site of stellar activity is seen towards our 
southernmost mapped region. At the base of the dense ridge 
seen in CS(l-O), we detect SiO emission in the v = 1,2,3 
modes of the (1-0) line. Here we are probably tracing the 
radiationally excited SiO emission from the variable star 
V5357Sgr. Spectra for the SiO emission are displayed in 
Figure [A2l 

We have also taken position switched deep pointings to- 
wards G5. 89-0. 39 B (DP-4) and towards the northern most 
point of the arm NE of G5.89-0.39 (DP-3). We have esti- 



mates of the mass and density from the CS isotopologue 
emission towards G5. 89-0. 39 B. We do not detect any iso- 
topologue emission from the extended arm position (DP-3) 
so we estimate the gas parameters under the assumption 
that the CS (1-0) emission is optically thin. The results of 
the compact-source analysis are displayed in Table [2] 

4.4 HESS J1800-240C / G5. 71-0.08 

The CS (1-0) emission detected towards HESS J1800-240 C 
shows two emitting regions along the line of sight, one with 
Vlsr ~ — 25kms _1 and another with Vlsr ~ 13kms _1 . 
This dense CS core is slightly offset from the peak of the 
TeV source and may provide a suitably dense region for 
CR interactions. The CS emission is also spatially matched 
with the SNR candidate G5. 71-0.08 and the 1720 MHz OH 
maser at Vlsr = 13kms _1 ([Hewitt fe Yusef-Zadehll2009h . 
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Figure 6. Integrated intensity maps and spectra for the southern HESS J1800-240B region. In all panels, black/white crosses x indicate 
the positions of G5. 89-0. 39 A (left) & B (right) and thick black contours are the H.E.S.S. TeV emission (4-6cr). White contours are used 
to highlight the molecular emission from the image. Minimum contour level for: CS (1-0) is 2.08Kkms _1 (~8 Trms); C 34 S (1-0) is 
0.32Kkms- 1 (-2T RMS ); 13 CS (1-0) is 0.24Kkm S - 1 (-2T RMS ); HC 3 N(5-4) is O.STKkms- 1 (-3T RMS ); CH 3 OH-I is O^Kkms" 1 
(~5Trms) an d SiO (1-0) is 0.36Kkms — 1 (~3Trms)- Surrounding the integrated images are spectra taken from both the mapping data 
and deep pointings, labelled according to their locations. 
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The OH maser may indicate that HESS J1800-240 C is 
tracing another SNR/molecular cloud interaction. Figure 
is an integrated intensity image of the HESS J 1800- 240 C 
region compar ing our CS ( 1 -0) em ission and the CO (2-1) 
emission from Fukui et al. (2008). Spectra from our C S 
(1-0) and CS(l-O) dee p pointings, the CO (2-1) s pectra 
from [Fukui et all (|2008l ) and CO (1-0) spectra from iLisztl 
(|2009l ) at the position of the 1720 MHz OH maser are also 
presented. 

The spectra presented in Figure [7] reveal line- wings in 
the CO (2-1) and (1-0) profiles, with Vlsr -0-8 and 20- 
30kms -1 , indicative of turbulent conditions. The CS (1- 
0) profile also shows some structure in these Vlsr ranges, 
which could be interpreted as multiple clumps along the line 
of sight, or evidence for line wings. There is also an absorp- 
tion dip at Vlsr ~ lOkrns -1 . Similarly, absorption features 
are present i n the 1720 MHz OH spectrum at Vlsr = 7 and 
-25kms _1 ([Hewitt fe Yusef-Zad eh 2009), and may also be 
evident in the CO emission. This is likely due to absorp- 
tion against a continuum source. Interestingly the width of 
the CS(l-O) line profile is narrow, with FWHM lkrns" 1 . 
This could suggest the central region of the core is rather 
quiescent. 

The two isotopologues of CS are detected in our deep 
position switched observation towards the new OH maser 
position and allow the mass and density to be estimated. 
For this CS core we estimate a mass of 190 M and a H2 
number density of 1.7xl0 5 cm -3 . Finally, we note that this 
mapped region associated with HESS J1800-240C has less 
exposure (a factor 1.5 to 3 lower) than the HESS J1801-233 
and HESS J1800-240 B regions. 



5 X-RAY & TEV GAMMA-RAY EMISSION 
AND DISCUSSION 

5.1 HESS J1801-233/NE Shocked Cloud 

The CS (1-0) emission seen towards the north-eastern 
cloud traces an extended dense region of molecular ma- 
terial. Using the emission within a region defined by 
/ T mb dv ^ lOKkms -1 (~ 40 T RM s) and finding the av- 
erage CS(l-O), C 34 S (1-0) and 13 CS (1-0) emission profiles 
within the region, we can estimate the gas parameters as an 
extended source. We apply this strict limit on the CS (1-0) 
emission to only accept pixels containing CS, C 34 S and 13 CS 
emission, and not bias the optical depth result. Under these 
conditions, we estimate the CS emission fills an ellipse on 
the sky with physical dimensions of 3.1 x 1.4 pc (at the esti- 
mated 2kpc distance). This elliptical region is indicated on 
FigureO Additionally, we assume the volume of the emission 
is defined by a prolate ellipsoid, where the third axis (in the 
z direction) is the same as the minor axis of the ellipse on 
the sky. In this way, the ellipsoid has radii 3.1 x 1.4 x 1.4 pc. 
Using the method outlined in S ection [3l we calculate the 
mass using iNicholas et al\ (|2010l Equation 3) and include 
an additional factor ?7mb/?7xb = 0.77 to acco unt for the Mo- 
pra e xtended beam efficiency ?7 x b = 0.43 ([Urquhart et all 
l201Ch . This gives an estimated extended source mass of 
5.6 x 10 4 M and a H 2 number density 6.1 x 10 4 cm" 3 . The 
extended source mass and density estimates are presented 
in Table [3 



We compare our CS gas mass to other mass esti- 
mates of the region. Over the 0-12kms _1 velocity range, 
the CO (1-0) mass has been estimated to be ~ 2 x 
1O 4 M , and ov er the wider 0-25 km s 1 velocity range 

- 5 x 10 4 M (lAharonian et all l2008bT ). From CO (3-2) 
data lArikawa et all (| 19991 ) estimated there to be 2 x 10 3 M 
of shocked gas in the NE cloud and 4 x 10 3 M© of unshocked 
gas. Additionally a conservative lower limit to the mass of 
gas traced by NH3 inversion transition s has been estimated 
to be > 1300 M ([Nicholas et qZ]|201Ch . 

At a glance the mass derived from CS is larger than 
previously published estimates. It is at least 10 x the lower 
limit obtained via NH3 observations. This is likely due to 
the considerable uncertainties in molecular abundance ra- 
tios, which are usually the largest sources of systematic error 
in mass/density calculations. While an order of magnitude 
mass discrepancy is not ideal, it has been shown for a sample 
of molecular cores, both hot and cold, that the core masses 
have equal or higher values from CS observations than those 
found via NH3 ; in some cases dis crepancies up to two orders 
of magnitude ([Zhou et all [l989). We suggest the CS mass 
be an upper limit as the abundance in a shocked source is 
likely to be lower than that we have assumed applies in an 
undisturbed, dense molecular cloud. The fact that the mor- 
phological extent of the CS emission is similar to that in CO 
however, suggests that a quite high fraction of the molecular 
gas is in a high density state. 

The X-ray morphology of W2 8 is described as cen- 
tre filled (|Rho fc Borkowskil [2002) whereby X-ray emis- 
sion lies towards the centre of the SNR. Interestingly, 
there is a brightening of the X-ray emission towards the 
N E of the SNR. This so called 'ear' feature, labelled 
bv iRho fc Borkowskil (2002), is found at the boundary of 
the pre- a nd post-shocked mol ecular cloud seen in CO (1-0) 
and (3-2) ([Arikawa et aZj|l999h . 

XMM-Newton X-ray observations towar ds the NE 
shock ed cloud have been discussed in lUeno et~all 
(2OO3]) f ollowing earlier observa tions with ASCA and 
ROSAT jRho fe Borkowskil |2003 ). Here, we have looked 
at additional XMM-Newton archival observations totalling 
^56 ks of observation time. Data were reprocessed using 
the XMM SAS software suite (with standard quality cuts) 
and major proton flares were filtered out. Observation IDs 
used were 013 5742401 (5. 5 ks po st filter), 0145970101 (20 ks 

- analysed bv lUeno et all (|2003l )). and 0145970401 (30 ks). 

The X-rays result mostly from th ermally heated ga s 
(post shocked) with kT = 0.3 keV (|Ueno et all I2003T ). 
Figure presents the XMM-Newton full EPIC camera 
(PN+MOS1+MOS2) image in the 0.2 to 10 keV band cor- 
rected for exposure in the blue channel. The red and green 
channels are the post-shocked CO (3-2 ) and pre-shocked 
CO (1-0) gas from 1 Arikawa efall (pL999). The right panel 
of Figure [8] is a zoom of the X-ray 'ear' feature, revealing 
the boundary between the X-ray emission and the shocked 
gas. Also seen at this boundary are many 1720 MHz OH 
masers (jClaussen et all Il997l groups OH E & OH F, la- 
belled on Figure 2}. Interestingly, adjacent to the OH maser 
groups E and F we also detected SiO (1-0) at the 2Trms 
level, further enhancing the evidence that there is shocked 
and disrupted dense gas in the region. 

The cluster of masers (OH E & OH F) are possibly 
associated with the X-ray 'ear 7 (Figure [8}. What is particu- 
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Figure 7. Left: Integrated intensity image of CS (1-0) emission towards HESS J1800-240C, which contains the recently detected 
1720 MHz OH maser (G5. 7-0.0) by Hewitt &; Yusef-Zadeh (2009). T he ma ser position is indicated by the cross +. Blue dashed contours 
are Nanten CO (2-1) integrated intensity contours from fFukui et all (|2008h . White contours are used to highlight CS (1-0) emission, only 
for the region of highest e xposure (indicated by the bla ck dashed b ox. Minimum white contour level is 0.88 K km s -1 (~2Trms)- Right: 
Spectra for the CO (2-1) ( Fukui et oZ]l2008h CO (1-0) (lLisztll2009h . CS (1-0) and C 34 S (1-0) for the maser position. Spectra for the CS 
isotopologue emission are taken from the Mopra deep pointing at the location of the OH maser. For the CS (1-0) emission, the vertical 
scale has been clipped (peak temperature 2.5 K) to show the evidence for line-wings and absorption in the profile. 



larly intriguing is that the OH E & F clusters of 1720 MHz 
OH masers appear to lie on both sides of the shocked gas 
(Figure [8] red channel) with respect to W28. This may re- 
sult from projection effects which are hard to disentangle or 
the possibility of an additional or reverse shock propagat- 
ing inward towards W28. The OH E & F maser clusters lie 
in the 10-15 km s _1 velocity range (Figure [4]). This suggests 
that the post-shocked X-ray emitting gas may lie behind 
most of the NE dense molecular cloud, which peaks in the 
5-10kms _1 velocity range (Figure SJ. 

The broad molecular line emission observed towards 
the NE cloud would likely result from non-thermal com- 
ponents. The additional kinetic energy required to produce 
the broad line emission can be calculated. We determine the 
non-thermal energy as VKkin = l/2M(Avki n ) 2 , where M is 
the mass of the broad line gas, and A^km is the line FWHM. 
Using the average CS (1-0) line FWHM=7.3 krns" 1 and the 
average extended mass upper limit M = 5.6 x 1O 4 M , 
we obtain Wkin < 3.8 x 10 49 erg, which is a few percent 
of the total 10 51 erg released in a typical SNR explosion. 
This upper limit is in agreement with our earlier estimate 
of the kinetic energy deposited into the dense NE cloud 
(Wcin > 7 x 10 47 erg) fr om our earlier NH3 observations. 
Similarly, lArikawa etal\ (|l999h estimated - 3 x 10 48 erg of 
kinetic energy has been deposited into the 2000 M traced 
by shocked CO (3-2). 

The impact of CRs and SN R shocks on molecular cl ouds 
has been discussed extensively. Ah aronian et al. dis- 



cuss a molecular cloud (MC) being overt aken by a SNR 
shock. The 'crushed cloud' model from lUchiv ama et all 
(2010) discusses the re-acceleration of particles within a 
molecular cloud , but i s mainly applicable to GeV energies. 
lYamazaki et all (|2006l ) discuss three cases for VHE emis- 
sion from old or evolving SNRs and molecular clouds (MC). 
This includes (i) emission from the SNR shock alone, (ii) 
the SNR shock colliding with a MC and (iii) CRs accel- 
erated a distance away from a MC and diffusing through 
the interstellar medium (ISM) to illuminate the cloud. 
This latter case of CR a cceleration and di ff usion is also 
more re cently discussed in Abd o et all (|2010h : iGabici et al. 
(|2010h : lOhira et all (l201l[). w ith discussion specifically to- 
ward W28. iFuiita et all (|2009h discuss the scenario of CRs 
accelerated by a SNR s hock nearby or ju st outside the W28 
MCs. Recent work by llnoue et all (|2010l ) models the prop- 
agation of middle aged SNR shocks into a cloudy ISM, re- 
vealing the presence of multiple secondary shocks with to- 
tal surface area comparable to the primary shock. These 
shocks are able to a ccelerate particles to a maximum en- 
ergy £ max > 10 TeV (jlnoue et q/.ll2010h . Overall, all models 
can account for the TeV and GeV emission seen towards 
W28 and other similarly aged SNRs. They clearly favour 
the hadronic emission mechanism for TeV emission. 

Under a hadronic emission mechanism, TeV photons 
result from CR collisions producing tt° mesons. The 7r°s 
then decay into the TeV photons. The leptonic emission 
mechanism produces TeV photons by inverse Compton pro- 
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Figure 8. Left; Three colour image for the HESS J1801-233 region. The red channel is JCMT CO (3-2) integrate d from -40 to 
40 kms" 1 ( Arik awa et all [19991 ). The green channel is the Nobeyama CO (1-0) emission integrated from 4 to 9 kms -1 ([Arikawa et all 
1999). The blue channel is the XMM-Newton 0.2 to lOkeV [phcm 2 s - 1 ] image (Gaussian smoothed). White dashed lines indicate the 
limits of the CO images. Thick white dashed contours are H.E.S.S. TeV emission (4 Sz 5cr levels). Yellow solid contours are Mop ra SiO (1-0) 
integrated from -5 to 20 kms -1 , as seen in Figure[3] White crosses + indicate the pos itions of 1720 MHz OH masers from Clausse n et all 
fll997f ). Right: Zoom of the X-ray ear and OH maser clusters OH E and OH F from lClaussen et all dl997h . 



cesses and/or electron non-thermal Bremsstrahlung. How- 
ever, these leptonic processes are suppressed in evolved 
SNRs due to electron cooling via synchrotron losses at much 
earlier epochs. 

TeV gamma-r ays are observed w ith energies in the 0.3- 
3 TeV range ([Aharonian et a/.ll2008bl ). In the hadronic sce- 
nario, typically 17 percent of the CR particle's energy is 
transferred to gamma-rays, therefore 0.3-3 TeV gamma-rays 
trace a parent population of ~ 2 — 20 TeV CRs at the source. 
The striking spatial overlap between the TeV emission and 
the dense gas towards HESS 1801-233, suggests that CRs 
are reaching the dense portions of the molecular cloud, 
and could help provid e constraints for CR diffusion mod- 
els (jGabici et aZj|2007h . However present H.E.S.S. observa- 
tions do not have sufficient angular resolution (~ 10 arcmin 
FWHM) to test whether the CR spectrum varies towards 
the interior of the dense cores. Future TeV instruments with 
sufficient angular resolution (arcmin) are required before 
such spectral features may become apparent and help to 
discriminate between models of re-acceleration and diffu- 
sion. A future TeV i nstrument, such as the Chere nkov Tele- 
scope Array (CTA) (|The CTA Consortium ll2010h . will have 
both improved angular resolution and sensitivity. The ex- 
pected angular resolution could reach a few arcmin FWHM 
and the energy sensitivity is expected to be a factor of 10 
better than H.E.S.S., achieving ^few x 10 -14 erg cm -2 s _1 
in rsj 50 hr. Together, these improvements will permit TeV 
emission maps and detailed spectra with comparable spatial 
scales to the molecular gas maps and provide a platform to 
test CR diffusion models and probe re-acceleration models. 

The expected flux of TeV emission from a molecu- 



lar cloud impacted by Galactic CRs (GCR) with particle 
flux following an integ ral power law E~ X ' Q can be found 
from Ah aronian 1 (|l99ll Equation 10): 

J F(^^ 7 ) = 2.85xl0- 13 ^Te 1 v (Jr 1 ) fcR [phcm-'s^Kl) 

where M5 is the mass of the cloud (in 10 5 M©), d^ pc is the 
distance to the cloud (kpc) and kcK is the CR density en- 
hancement factor above the l ocal solar- system value. T he ob- 
served TeV spectral indices ([Aharonian et aZ.ll20Q8bh imply 
an incident proton particle spectral integral index of F ^ 1.6 
(in a power law particle spectrum oc E~ r ) similar to GCRs. 
Thus Equation [1] should be sufficient for initial estimates of 
the TeV fluxes arising from parts of the clouds. We note that 
further discussion of other systematic uncertainties concern- 
ing the predicted TeV fluxes from cloud clumps is left for 
section [6] 

The CR enhancement factors have b een estimated to- 
wards three of the four TeV sources in ([Aharonian et all 
2008b, Table 2). With this in mind, we can estimate the 
gamma-ray flux for emission from the dense cloud compo- 
nents traced by CS (1-0) emission. Using the NE shocked 
cloud extended source mass estimate from Table [3] we es- 
timate F(> ITeV) = 5.1 x 10~ 13 phcm -2 s" 1 . The flux 
results for this, and other CS regions, are summarised in 
Table [H We note that this flux is ~ 15% higher than that 
currently detected by H.E.S.S. This further supports our 
case for treating the CS mass estimate as an upper limit 
given the uncertainties in the CS abundance. The general 
picture though is that the predicted TeV flux is detectable 
and can be well studied by CTA in a reasonable observation 
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time (10-50 hr). CTA could then discriminate between the 
TeV emission from within the dense cloud component vs. 
that of the entire cloud as seen by H.E.S.S.. 

5.2 HESS 1800-240 A & B 

The southern TeV sources HESS J1800-204 A and B are also 
associated with dense gas which could provide a suitable 
target for CR interactions. The dense regions of gas spa- 
tially consistent with HESS 1800-240 A are treated as el- 
liptical extended sources filling prolate ellipsoids with radii 
1.4 x 0.7 x 0.7 pc for the gas nearby G6. 225-0569 and 
2.8 x 1.4 x 1.4 pc for the gas nearby G6. 1-0.6. Similarly, the 
CS gas seen towards HESS J1800-240 B is treated as three 
separate extended sources. The large clump associated with 
the G5. 89-0. 39 Hn regions is assumed to fill a prolate ellip- 
soid with radii 2.4 x 2.1 x 2.1 pc. The CS emission from the 
NW arm, extending from V5561 Sgr, is assumed to fill an 
ellipsoid with radii 3.1 x 1.1 x 1.1 pc, and the SE molecular 
ridge is assumed to fill an ellipsoid with radii 3. 5 x 1.4 x 1.4 pc. 
These elliptical regions are indicated on Figure [2] We detect 
only CS (1-0) emission from the NW arm and SE ridge re- 
gions, therefore, we assume the emission to be optically thin. 
The individual estimates for the mass and density of all the 
extended regions towards HESS J 1800- 240 A & B are sum- 
marised in Table [3] 

The three extended regions towards HESS J1800-240 B 
suggest that 7.1 x 10 4 M of material is contained in 
dense clumps. The HESS J1800-240 B southern cloud is 
estim ated to have 10 5 M fi b ased on the CO (1-0) emis- 
sion (jAharonian et a l 2008b). Based on our extended CS 
clump estimates, we find that ~ 70% of the southern clouds 
mass is contained in the dense clumps surrounding the 
G5.89-0.39 Hn regions. 

As these three dense extended regions all lie within the 
HESS J 1800- 240 B source, the improved angular resolution 
of CTA could potentially resolve this emission into three sep- 
arate sources. In this case, the expected gamma-ray flux can 
be found via Equation [1] The expected TeV flux for photon 
energy E > 1 TeV for each of the dense cloud components 
are summarised in Table |4] Typically the expected flux is 
10~ 13 to 10~ 14 ph cm s , which could be detectable and 
well studied by CTA in a reasonable observation time (10- 
50 hr). As for the NE cloud, we find that the total expected 
flux from the three HESS J1800-240 B sources is higher than 
the detected level of TeV emission (by a factor of ^2). Again, 
this could be explained by our choice of CS abundance being 
overestimated in this region which may be expected given 
likely internally shocked nature of the G5. 89-0. 39 Hll com- 
plex. Additionally, our cloud- averaged CR density ken may 
not reflect smaller scale variations. 

The question of the origin of the CR population towards 
HESS J1800-240 A and B is still uncertain, although diffu- 
sion of CRs from W28 seems to reasonably wel l explain the 
observ e d GeV to TeV gamma-ray spectra (e.g. iFuiita et al. 
(2009h; lAbdo et all (|2Q10h : iGabici et all (|2010h : iLi Chen 
(2010)). However, additional CRs accelerated inside the lo- 
cal star formation regions (from G5. 89-0. 39 in particular), 
and another potential SNR towards HESS J18 00-240 C (dis- 
cussed shortly) remains a possibility (see e.g. lAraudo et al\ 
(2007)). Our observation of the arcmin-scale dense NW arm 
and SE ridge features in the molecular cloud may signal 



additional sites of disruption which could arise from in- 
ternal (star formation) and/or external (e.g. SNR shocks) 
forces. Further tests of these scenarios can come from future 
arcmin-scale studies of gamma-ray spectra throughout the 
clouds. 



5.3 HESS J1800-240C 

Towards HESS J 1800- 240 C, we have discovered a dense 
molecular core. Unlike the other TeV sources in the W28 
region, the dense gas lies slightly offset from the peak of 
the TeV emission. The CS gas is found in a dense, com- 
pact core, which has a narrow line width, FWHM lkrns" 1 . 
Additionally, there is absorption in the CS line profile at 
~10kms _1 , which may also be evident in the CO (1-0) and 
(2-1) line profiles (see spectra in Figure [7]). 

Bas ed on the CO (1-0) em ission from the NANTEN 
survey (Mizu no & Fukui 2004), we estimated the mass 
of the CO cloud. We assumed a CO abundance of Xco 
— 1.5 xlO 20 and assumed a spherical emission region 
with radius 0.1° centred on the HESS J1800-204C TeV 
source. We included emission from the 0-20kms _1 range, 
whic h was selected to be c onsistent with the ranges used 
in Aharonian et al. | (|2008bh . Under these assumptions, we 
estimated 1.4xlO 4 M of CO gas towards HESS J1800- 
240 C. We estimated the CR enhancement factor ken = 35, 
using Equation [T] based on the CO mass a nd the detected 
TeV flux following lAharonian et all (|2QQ8bh . Assuming the 
value of kcR = 35 we determined the expected TeV flux 
from just the dense CS core. The result is shown in Tabled 

The determined CR enhancement of 35 is higher than 
towards the other TeV sources. Enhancements of this level 
are expected in the regions of CR accelerators, which could 
be expected as this region could be associated with both 
the W28 SNR as well as the local SNR candidate G5.71- 
0.08 which may be producing the OH maser. Under a purely 
hadronic emission scenario, the expected TeV flux from just 
the dense CS core is F(E ^ 1 TeV) - 2 x 10" 15 ph cm" 2 s -1 . 
This flux level would likely push the reasonable detection 
limits for CTA. 



6 SUMMARY AND CONCLUSIONS 

Using the Mopra 22 m telescope, we have conducted 7 mm 
molecular line mapping covering dense and disrupted gas 
toward the W28 SNR field. We have followed up our previ- 
ous 12 mm line study of the molecular clouds towards the 
TeV gamma-ray peaks observed by the H.E.S.S. telescope, 
with higher resolution targeted observations towards the NE 
cloud/shock interaction, HESS J 180 1-233 and the bright and 
energetic UC-Hn regions in the southern cloud HESS J 1800- 
240 B. 

Sites of shocks, outflows and disruption are revealed 
with SiO. In the majority of cases, SiO detections are to- 
wards sites of stellar activity and Hn regions, however, SiO is 
also detected in the NE cloud/SNR shock interaction region 
with broad line profiles. Interestingly, in this shock region 
the SiO emission is bounded by clusters of 1720 OH masers, 
indicating the downstream direction of the shock is towards 
W28. Other lines detected are 44 GHz CH 3 OH masers and 
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Table 4. Summary of the CR enhancement factors, /ccr, and 
the predicted TeV fluxes from the extended dense regions listed 
in Table [3] The predicted fluxes assume a CR enhancement factor 
determined in lAharonian et all (120081)1 ). except for HESS J1800- 
240 C, kcR has been calculated here. We treat the flux from the 
NE shocked cloud as an upper limit, as the estimated mass is an 
upper limit. 

Region k C R Expected Flux F(E ^ 1 TeV) 

[phcm -2 s _1 ] 



NE Shocked Cloud 


13 


< 5.1x10 


-13 


G5.89-0.39 Hn Region 


18 


7.0x10- 


13 


G5.89-0.39 NW Arm 


18 


5.4x10- 


14 


G5.89-0.39 SE Ridge 


18 


1.6x10- 


13 


G6. 1-0.6 Region 


14* 


6.5x10- 


14 


G6.225-0.569 Region 


14* 


1.6x10- 


14 


HESS J1800-240C 


35t 


1.9x10- 


15 



* kcR from I Aharonian et all (|2008rJ ) has been re-scaled 
assuming a d = 2 kpc. 

^ kcB, calculated here us ing; the CO (1-0) emission from 
iMizuno fe Fukuil (120041 ) 

■f- Assuming the compact-source size mass indicated in Table [2] 



cyanopolyynes HC n N (n = 3,5,7) which are again prefer- 
entially associated with star formation and Hll regions. 

The CS cores are typically found towards CO peaks 
and trace conditions ranging from quiescent cloud cores to 
Hn regions. The obvious exception is the NE cloud/SNR 
shock interaction region, which exhibits broad line emission 
FWHM > lOkms -1 in all detected lines. 

Based on CS isotopologue ratios from our position- 
switched deep pointing observations, we are able to estimate 
the upper state CS column density. Temperature estimates 
from our previous NH3 observations allow the LTE molecu- 
lar hydrogen mass and density to be estimated. Our results 
presented in Table [2] assume a compact-source core radius 
r = 0.2 pc, however, scaling factors for various core radii are 
also included. 

We have estimated the mass of the extended dense NE 
cloud and several of the dense clumps in the southern 
clouds, assuming the CS isotopologue emission is contained 
within a prolate ellipsoid. Under these assumptions, we es- 
timate there to be 5.6xl0 4 M©, with H2 number density 
6.1 xlO 4 cm -3 (as an upper limit) contained towards the 
NE HESS J1801-233 TeV source. We find that - 4 x 10 49 erg 
of kinetic energy is required to produce the broad emission 
in the CS line profile towards the NE cloud, which is a few 
percent of the typical 10 51 erg of kinetic energy released in 
aSNR. 

For the southern cloud towards HESS J1800-240B we 
estimate that ~ 70% of the gas mass is contained in dense 
clumps surrounding the Hn regions. The southern clouds 
show no signs of external disruption, providing no evidence 
that the W28 SNR shock has influenced the region. 

Although our CS masses are treated as upper limits, we 
note that a somewhat independent estimate of the CS mass 
for clumps not part of the NE Shocked Cloud comes from the 
virial theorem. For these cases the fact that the virial mass 
is similar to the LTE mass within a factor 2 to 3 suggests our 
choice of CS abundance is adequate for order of magnitude 
estimates. For the NE Shocked Cloud the virial theorem will 



likely not apply. However, as discussed in section [5TT1 the NE 
Shocked Cloud mass derived from our CS observations is 
10% larger than the mass from CO observations, indicating 
that our CS abundance may be slightly overestimated here. 

Additional systematic uncertainties arise from the un- 
known spectrum of CRs entering the dense clouds, and, pen- 
etrating their interiors as a result of energy-dependent diffu- 
sion. Equation [T] assumes that a CR spectral integral index 
of r = —1.6 as expected for Galactic CRs (GCRs) impact- 
ing passive clouds some distance away. The spectral indices 
of the TeV emission towards W28 are in fact found to be 
similar to GCRs so Equation [1] should be adequate in this 
regard. 

The striking spatial match between the TeV emission 
and the dense molecular gas suggests that the CRs are able 
to penetrate some distance into the dense cloud core. How- 
ever, the energy depen dent diffusion of CRs into the cloud 
cores ([Gabici et a/.ll2007h could easily suppress the predicted 
emission for low energies E < 1 TeV below that suggested 
by Equation [TJ The level of suppression is strongly depen- 
dent on the largely unknown magnetic field structure and 
turbulence in the clump and is a topic beyond the scope of 
this paper. 

Overall, we would like to emphasise that the uncertain- 
ties discussed above clearly motivate the need for new 1 
arcmin angular resolution TeV observations (for example by 
CTA) in order to begin probing the diffusion properties of 
CR into dense molecular clouds. 

This work is part of our ongoing study into the molec- 
ular gas towards the W28 region, and follows our earlier 
12 mm line mapping campaign. Additional observations of 
the NE cloud SNR interaction region have been conducted 
with Mopra at 12 mm in order to extract detailed NH3 in- 
version spectra across the cloud core. Results of these ob- 
servations will be the focus of a future paper to further un- 
derstand the molecular gas properties in the interesting NE 
SNR shock/molecular cloud interaction region. 
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7mm line survey: W 28 field TeV sources. 19 
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Figure Al. Spectra for the line emission detected towards the 
IRAS 17589-2312 star formation core from mapping data. 

APPENDIX A: ADDITIONAL FIGURES AND 
TABLES 
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Figure A2. Spectra of the SiO (1-0) v = 1, 2, 3 lines coming from 
the V5357 Sgr Variable Star of Mira Cet type. 




2^s 18h02m00s 36s 01 ml 2 s 00 m 4 

Right Ascension (J2000) 



Figure A3. Image showing the relationship between the diffuse 
and dense gas in the NE cloud. The colour scale is the NANTEN2 
CO (2-1) -10 to 25kms _1 integrated intensity image. The black 
dased contours are the H.E.S.S. TeV emission (4 & 5 a levels). The 
white dashed contours are the CS (1-0) -5 to 20kms _1 integrated 
emission as seen on Figure[3] The solid cyan contours are the NH3 
(3,3) -5 to 2 0kms~ 1 integrated em ission (1,3,5,7,9,11,13 K km s _1 
levels) from lNicholas et al\ (l2010h . 



© 0000 RAS, MNRAS 000, 000-000 



20 B. Nicholas et al. 

Table Al. Line parameters from a Gaussian fit to the spectra, for lines detected in deep pointing observations. 



Region Peak T mb J T mb dv V LSR FWHM 

[K] [Kkms- 1 ] [kms- 1 ] [kms" 1 ] 



- DP-1 - 



CH30H 


0. 


,03 


± 


0. 


.003 


0.55 ± 0.09 


11.70 ± 1.072 


18.74 zb 2.29 


SiO (1-0) 


0. 


.09 


zb 


0. 


,003 


2.15 ± 0.12 


13.12 ± 0.395 


22.71 zb 0.98 


HCc>N (5-4) 


0. 


,04 


± 


0. 


,006 


0.61 zb 0.14 


12.18 ± 1.215 


15.23 zb 2.65 


13 CS (1-0) 


0. 


,02 


zb 


0. 


,003 


0.58 ± 0.13 


13.68 db 1.778 


25.66 zb 4.63 


C 34 S (1-0) 


0. 


,07 


zb 


0. 


,006 


1.39 ± 0.19 


12.30 ± 0.603 


17.80 zb 2.05 


CS (1-0) 


1 


.06 


± 


0. 


005 


20 46 + 16 


13 05 + 043 


18 09 + 11 














- DP-2 - 






CH30H 





,06 


zb 


0. 


,004 


0.98 ± 0.09 


5.67 dz 0.458 


15.31 zb 1.08 


SiO (1-0) 


0. 


.13 


zb 


0. 


,004 


2.29 ± 0.11 


6.10 ± 0.231 


15.99 zb 0.58 


HC 3 N (5-4) 


0. 


,10 


zb 


0. 


.007 


1.16 ± 0.12 


6.87 ± 0.361 


10.40 zb 0.84 


13 CS (1-0) 


0. 


,07 


± 


0. 


,006 


0.79 db 0.10 


8.43 db 0.352 


10.49 zb 1.04 


C 34 S (1-0) 


0. 


,20 


± 


0. 


,007 


1.99 ± 0.12 


7.55 ± 0.168 


9.60 zb 0.45 


CS (1-0) 


2. 


.37 


zb 


0. 


,007 


28.27 ± 0.13 


7.37 ± 0.015 


11.20 zb 0.04 














- DP-3 - 






CS (1-0) 


1 


,39 


zb 


0. 


,064 


3.28 ± 0.22 


8.98 dz 0.043 


2.21 zb 0.11 














- DP-4 - 








3 


.65 


-j- 


0. 


014 


4 Aa _l n no 


q 00 -\- 009 


1 1 c _L n m 

±.±0 _l_ \J .\J -L 


SiO (1-0) 


0. 


,25 


± 


0. 


,006 


8.98 ± 0.05 


12.35 zb 0.011 


19.56 zb 0.02 


HC 5 N (17-16) 


0. 


,59 


zb 


0. 


,005 


5.15 ± 0.17 


9.17 zb 0.152 


3.29 zb 0.53 


HC 3 N (5-4) 


4 


.12 


± 


0. 


,009 


2.07 ± 0.04 


9.09 zb 0.024 


3.56 zb 0.05 


13 CS (1-0) 


0. 


.40 


zb 


0. 


,013 


15.63 ± 0.07 


9.35 zb 0.005 


3.64 zb 0.01 


HC 5 N (16-15) 


0. 


.52 


zb 


0. 


,009 


1.53 ± 0.06 


8.97 zb 0.037 


3.10 zb 0.10 


C 34 S (1-0) 


1. 


.15 


zb 


0. 


,013 


1.72 ± 0.06 


9.29 zb 0.037 


3.67 zb 0.08 


CS (1-0) 


9. 


.50 


zb 


0. 


,015 


4.50 ± 0.09 


9.26 zb 0.022 


4.01 zb 0.06 














- DP-5 - 






C 34 S (1-0) 


0, 


,29 


zb 


0. 


,022 


0.31 ± 0.04 


12.97 zb 0.037 


1.01 zb 0.09 


CS (1-0) 


2 


.80 


=b 


0. 


,022 


3.34 ± 0.04 


13.04 zb 0.004 


1.12 zb 0.01 
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